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We describe a general approach to quark flavor tagging 
in polarized hadronic processes, with particular emphasis on 
semi- inclusive deep inelastic scattering. A formalism is intro- 
duced that allows one to relate chosen quark flavor polariza- 
tions to an arbitrary combination of final-state hadron spin 
asymmetries. Within the context of the presented formalism, 
we quantify the sensitivity of various semi-inclusive hadron 
asymmetries to the light quark flavors. We show that unpo- 
larized A's may allow one to measure strange quark and anti- 
quark polarizations independently. We also highlight several 
applications of our formalism, particularly to measurements 
intended to probe further the spin structure of the nucleon. 

PACS number(s): 13.88.+e,13.60.-r,13.85.Ni,13.87.Fh 



I. INTRODUCTION 

The interpretation of hard hadronic processes in terms 
of partonic (quark and gluon) degrees of freedom forms a 
vital component of modern research in particle physics. 
In particular, the tagging of parton quantum numbers, 
such as flavor and charge, using appropriately chosen fi- 
nal states allows important tests of both QCD and elec- 
troweak dynamics. In the remainder of this paper, we 
will focus primarily on the leptoproduction of hadrons in 
deep inelastic scattering (DIS), IN — > I'hX, where both 
initial states are polarized. Such measurements continue 
to provide important information on the spin structure of 
the nucleon |Q. However, the conceptual foundations of 
our approach and the formalism that we introduce per- 
tain to any hard hadronic process, and other applications 



are mentioned in Sec. IV 



The spin structure of the nucleon has received much 
attention in the past decade Q. The experimental fo- 
cus until rather recently has been on the leading twist 
structure function gi(x, Q 2 ), which is roughly propor- 
tional to the inclusive spin asymmetry on a longitudinally 
polarized target. However, processes where at least one 
hadron is detected in the final state offer several distinct 
advantages over the inclusive process alone §-|. In par- 
ticular, semi-inclusive reactions provide a direct probe 



of the flavor dependence of quark observables, allowing 
more stringent tests of hadron structure. Perhaps more 
importantly, semi-inclusive processes allow one to sepa- 
rate [||J5) contributions with definite charge conjugation 
symmetry. This is particularly interesting since charge- 
odd quantities are free of the axial anomaly. It has been 
recently suggested || that a comparison of the charge 
conjugation even and odd parts of parton helicity and 
transversity distributions may yield new information on 
the polarized gluon contribution to the nucleon spin. 

Given this importance of semi-inclusive measurements 
in hard hadronic processes, it is useful to have a for- 
malism that is independent of model assumptions. In 
Sec. Ill, we introduce a general method to extract cho- 



sen polarized quark distributions from an arbitrary com- 
bination of final state asymmetries. This new approach 
is particularly useful in light of the latest generation of 
existing [HERMES, Spin Muon Collaboration (SMC), 
TJNAF] and forthcoming [COMPASS, ELFE, HERA-p, 
BNL Relativistic Heavy Ion Collider (RHIC) Spin] ex- 
perimental efforts, which allow identification of a large 
number of final state hadrons. Within the framework of 
the presented formalism, we quantify the sensitivity, as 
a function of kinematics, of various hadrons to the light 
quark flavors. 



II. GENERAL CONSIDERATIONS 

Traditionally, quark flavor tagging has been applied 
almost exclusively to the current fragmentation region 
in deep inelastic processes. Since a separation between 
the current and target regions is usually regarded as a 
necessary criterion for such measurements, one imposes 
kinematic cuts in order to try to exclude the target re- 
gion. Our point is that such a separation is not pos- 
sible, even in principle, within a coherent and complete 
description of hadronic final states]] In perturbation the- 
ory, this follows from the existence of collinear singulari- 
ties with respect to initial state partons that generate a 
contribution to the target fragmentation region 0. We 
further argue that, generally speaking, there is no a pri- 
ori reason to minimize the effects of target fragmenta- 
tion. Indeed, it has been shown M that, independent of 
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1 Sophisticated hadronization models, such as the string and 
cluster models, reflect this fact. 
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a transverse momentum cutoff, all collinear singularities 
for sufficiently hard hadrons can be absorbed into so- 
called fracture functions || M^ n (x,Xf,Q 2 ), which give 
the joint probability of "finding" within the target N a 
parton i and hadron h. Measurements sensitive to target 
fragmentation offer complementary information on the 
nucleon state |^| and additional insights into a unified 
view [fiof of hadronic reactions. 

Quark flavor tagging is based on the idea of local 
parton-hadron duality (LPHD) |Tjj| , where it is assumed 
that the flow of quantum numbers at the hadron level 
tends to follow the flow established at the parton level. 
At leading order in perturbation theory (for infrared safe 
or factorizable parton quantities), this has the practical 
consequence that the parton-hadron correspondence is 
essentially one-to-one. Since quarks and gluons are not 
asymptotic states (due to confinement), it would be a 
mistake to conclude that we can ever measure, in effect, 
a primary parton, instead of merely an event property 
that is correlated with the primary parton. Within a 
complete description of hadronic processes, valid beyond 
the leading order, the relevant issue is the degree of corre- 
lation between the final state topology and the quantum 
numbers of the primary partons. 

It is therefore instructive to discuss briefly the nature 
of parton-hadron correlations. Inclusive hadron distri- 
butions in all hadronic processes, regardless of whether 
the collisions are hard or soft, are characterized by pro- 
jectile fragmentation regions separated by a central re- 
gion Jl^] . The key feature of the central region is that it 
is essentially independent of both projectiles and hence 
universal for all hadronic processes (at fixed invariant 
mass W). This has been confirmed in hadron-hadron 
collisions fL3|| , and indeed, hadron spectra in the central 
region are found to be very similar [ 14|] in pp, pp, photo- 
production, and low- a; DIS processes. These generic fea- 
tures of hadronic final states have been recognized [ ^5|Jl2| 
a long time ago to be a consequence of Lorentz invariance 
and short range correlations in rapidity. 

Using an approach based on short range rapidity corre- 
lations and LPHD, we have calculated Q the correlations 
of light hadrons with respect to the current quark and 
target remnant. Within our discussion, it is important 
to distinguish between forward (xf > 0) and backward 
{xp < 0) regions, which are defined strictly by kinemat- 
ics, from the current, target, and central fragmentation 
regions, which are never distinct, but represent varying 
degrees of correlation with the quark and remnant. As 
the correlations depend strongly on rapidity differences, 
our definition of the fragmentation regions coincides with 
their classification, within perturbative QCD, in terms of 
collinear and soft divergences. This suggests that it is fa- 
vorable to use hadronic variables such as xp or rapidity 
over the usual energy fraction z = E] l jE 1 , since the lat- 
ter variable cannot distinguish between the central and 
target regions (both dominate at low z). 

In our calculations, we find that at small \xf \ < 
0.1 — 0.2 and fixed Q 2 , both the current and target corre- 



lations decrease as W increases.]] This implies that most 
of the increase in hadron production with increasing W 
occurs in the central region, a result well known from 
hadronic phenomenology. Indeed, in DIS at constant Q 2 , 
the charged pion fragmentation functions D 7T (xf, W) rise 
linearly with hiW at small |a^| < 0.1, but remain con- 
stant at larger values of \xf\, both in the forward and 
backward regions |lS|| . This is an important result, be- 
cause it illustrates that independent of the target rem- 
nant, there are kinematical effects in the forward region 
that cannot be described by fragmentation functions de- 
pending on z and Q 2 alone, as is usually assumed. For 
increasing values of W , while the current and target re- 
gions become better separated kinematically, the fraction 
of hadrons that are strongly correlated with the current 
quark decreases (roughly as 1/lnVF). Even in the cen- 
tral region, however, the correlation between an arbitrary 
final state and the flavor of the current quark remains sig- 
nificant due to the finite number of accessible flavors. In 
the next section, we introduce a new formalism that al- 
lows one to quantify the above correlations and to relate 
chosen hadron-level asymmetries to the polarization of 
partons involved in the hard scattering process. 



III. FORMALISM 

Hadronic final states in deep inelastic processes are 
typically analyzed in terms of the usual fragmentation 
functions D|* H , which parameterize the probability den- 
sity for producing hadrons h, given an initial state parton 
of flavor i (i = u, u, d,d, . . .). The corresponding formu- 
las for polarized and unpolarized scattering, at leading 
twist, are well known. However, fragmentation functions 
are not the most useful quantities in the context of quark 
flavor tagging. Instead, we define the quark flavor purity 
Pj 1 as the probability that a quark of flavor i was probed 
by the virtual photon, given a final state h. In terms of 
the usual quantities, we can write 



ph 



D 



f e i 



■v 



(1) 



where are the quark charges, qi — qi{x,Q 2 ) are helic- 
ity averaged quark distributions, Df = D^xftQ 2 , x ) are 
generalized fragmentation functions]^] and Yli=f j = 
1 holds for any final state. The purities not only have the 



Several quantitative examples were given in Ref. @ on 
the basis of scaling violations [^| in inclusive charged parti- 
cle momentum spectra in e + e~ collisions. However, as dis- 
cussed in Refs. |fl7yi8| , these scaling violations are actually 
due to threshold production of charmed particles, where, as 
expected, the region of scaling extends down to lower values 
of xf as W is increased. 
3 For convenience, we integrate over the hadron's pr and (f>. 
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physically meaningful interpretation mentioned above, 
but also are the relevant quantities for performing a po- 
larized quark flavor decomposition from semi-inclusive 
asymmetries. 

As a specific example, we consider semi-inclusive deep 
inelastic scattering with both beam and target longitudi- 
nally polarized, where the photon is the exchanged boson. 
In this case, for hadrons h in the final state, one measures 
the virtual photon-nucleon asymmetry 



A h {x,x F ,Q 2 ) = 



'1/2 



'3/2 



(2) 



'1/2 



'3/2 



where <rjy 2 (3/2) ( x i x f,Q 2 ) denotes the cross section when 
the projection of the total angular momentum of the 
photon-nucleon system is 1/2(3/2). We can then write 
any such asymmetry (at leading twist) in terms of only 
purities and quark flavor polarizations (Aq/q)f. 



A" 



i=f,f 



(3) 



Following the arguments in the previous section, Eq. (|^) 
holds in both the forward and backward regions, for suf- 
ficiently hard hadrons. The derivation of Eq. (|^) depends 
crucially on the independence of the fragmentation pro- 
cess on the quark polarization. Although this is usually 
stated as an additional assumption, it is actually a conse- 
quence of the parity invariance of the strong interaction: 
for inclusive production of unpolarizcd hadrons, there is 
no pseudovector observable in the final state that could 
in principle couple to the quark polarization. It follows 
that, for example, the presence of target fragmentation 
does not dilute the asymmetry. In general, there is no 
reason to minimize the contribution of target fragmen- 
tation, as long as the correct purities are used in the 
analysis. 

The form of Eq. (||) allows us to combine several inde- 
pendent asymmetries into a column vector, yielding the 
matrix equation: 



q 



(4) 



When the matrix elements of V are determined (from 
data or simulations), Eq. (^) contains all information 
needed to extract specified polarized quark distributions 
from a chosen set of semi-inclusive asymmetries. The in- 
clusive asymmetry, where only the scattered lepton is de- 
tected, also contains useful information and can be triv- 
ially incorporated within the purity formalism: 



p-inci 



=/,/ 



1 el qv 



(5) 



Within the presented approach, then, the extraction of 
polarized quark flavor distributions is reduced to opti- 
mization, via choice of final (and initial) state hadrons 



and kinematics, of the purity matrix. Several methods 
of optimization are suggested in Ref . ||] . 

We have studied the kinematical dependence of puri- 
ties for light quark flavors and hadrons, on proton and 
deuteron targets, using lepto-6.1 pQ] for event gen- 
eration and jetset-7.4 J2TJ for hadronization. Since 
the purities are independent of the quark helicities, an 
unpolarized simulation can be used to generate them. 
We have used the HERMES kinematics and a paramet- 
ric model ^2) of its acceptance in our simulations (see 
Rcf. for details), applying standard DIS selection cuts 
of 0T> 1 GeV 2 , W 2 > 4 GeV 2 , and y < 0.85. The 
following numerical results therefore apply primarily to 
the forward region, as HERMES detects mostly forward 
hadrons. As expected from the form of Eq. ([[]), the puri- 
ties are relatively insensitive to the choice of hadroniza- 
tion model parameters and unpolarized parton distribu- 
tions. 
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FIG. 1. Purities Pf{x, \xf\ > 0.1) for typical final 
states with respect to the light quark flavors: it (solid), 
u (short-dashed), d (long-dashed), s (dotted), and s 
(dot-dashed). The d-quark contribution never exceeds 10% 
and is not shown. 

In Fig. [l], we plot quark flavor purities for typical fi- 
nal states as a function of the Bjorken-x variable, us- 
ing a deuteron target (results on the proton are qual- 
itatively similar). All quark flavor labels refer to the 
proton, using isospin symmetry. The cut \xf\ > 0.1 
is applied to suppress very central hadrons. Using the 
purities shown, it is straightforward to quantify the sen- 
sitivity of semi-inclusive asymmetries to the light quark 
and antiquark flavors. We summarize our general conclu- 
sions here. As expected from charge counting, positively 
charged hadrons are primarily sensitive to the u-quark 
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[see Figs. |T](a),(c)]. Given that the inclusive asymmetry 
on a proton target essentially probes the combination 
Au + Au, one can reasonably constrain the u- and u- 
quark polarizations from a combined analysis of inclusive 
and positively charged hadron data alone. Since PJ is 
dominated by it-, u-, and <i-quarks at small x < 0.1 [see 
Fig. 0(b)], including A* in the above analysis allows 
one to extract Au/u, Au/u, and Ad/d in the sea region. 
In addition, as only u- and d-quarks contribute signifi- 
cantly to PJ in the valence region (x > 0.2), one can 
extract Au/u and Ad/d (at large x) using only A* and 
A 17 . We conclude that a combination of asymmetries 
using copiously produced hadrons can be used to extract 
Au/u, Au/u, and Ad/d. On the other hand, it will be a 
challenge to measure the polarization of <i-quarks, since 
the corresponding purities never exceed 10% for the final 
states shown in Fig. [l| We suggest combining Ad with a 
measurement of Ad v (= Ad — Ad) using the pion charge 
difference asymmetry method of Ref . H .[] 
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FIG. 2. Purities Pj(xf,x < 0.1) for strange hadrons with 
respect to the light quark flavors: u (solid), u (short-dashed), 
d (long-dashed), s (dotted), and s (dot-dashed). The d-quark 
contribution never exceeds 10% and is not shown. In (b), the 
xf variable refers to the daughter K + . 

We have also studied the sensitivity of semi-inclusive 
asymmetries to polarized strange quarks and antiquarks 



4 The derivation in Ref. j| assumes charge and isospin con- 
jugation symmetries for the fragmentation functions, which 
are formally violated in a complete QCD description of DIS. 
Phenomenological consequences are studied in Ref. p3|. 



in the nucleon. In Fig. we plot the purities of selected 
strange hadrons as a function of xp in the sea region 
(x < 0.1). The K~ is often regarded as an interesting 
probe of the nucleon, as it is an "all sea" object. We 
indeed find a large sensitivity to sea quarks. As expected 
from charge counting, the u-quark dominates over the s- 
quark contribution, independently of xf [see Fig. ^(a)]. 
On the other hand, selecting charged kaons from <f> meson 
decay provides improved sensitivity to s- and s-quarks 
[see Fig. |](b)]. For example, using the cuts x < 0.1 and 
x F > 0.1, we find that Pf + Pf ~ 40-45%. The meson 
may therefore be a useful probe of the polarized strange 
sea, provided one applies kinematical cuts to suppress 
elastic production. 

Semi-inclusive A (A) asymmetries arc predominantly 
sensitive to u- {u- and u-) quark polarizations, except at 
rather large values of [see Figs. |](c),(d)]. However, 
we find that the A is sensitive, at the 20% level, to s- 
quarks at slightly negative xf < —0.2, where the event 
rate is still reasonably large. The key point is that the 
A asymmetry in this region is determined essentially by 
Au/u, Ad/d, and As/s (and the corresponding purities). 
Since Au/u and Ad/d will be well constrained by charged 
hadron and inclusive asymmetries (as discussed above), 
A production in the backward region may be used to 
constrain the polarized s-quark sea. In this case, to a 
good approximation, one can write 



As 
s 



1 

w 



A' 



E 

i—u,d 



(6) 



Likewise, as the A event rate distribution in xf is rather 
broad, one obtains appreciable sensitivity to polarized s- 
quarks using only a modest cut on xf- In particular, for 
x F > 0.25 (and x < 0.1), we find that P^ ~ 55% and 
P A ~ pA _ 20%. Hence, given the contribution of u- and 
d-quarks, the A asymmetry in this region probes As/s. 
We therefore emphasize the importance of measuring A 
asymmetries in both the forward and backward regions. 



IV. OTHER APPLICATIONS 

Thus far, we have introduced the quark flavor purities 
with particular emphasis on their key role in interpret- 
ing electroproduction of unpolarized hadrons, where both 
beam and target are longitudinally polarized. However, 
it is clear that our considerations generalize directly to all 
sufficiently hard hadronic processes. In general, the puri- 
ties allow us to quantify, in a well defined way, the sensi- 
tivity of chosen hadrons to the flavor of quarks involved in 
the hard scattering process. For the special case of polar- 
ization in the initial state, the purities relate, within the 
framework of a compact formalism, hadron-level asym- 
metries to the polarization of initial state partons. For 
example, one may imagine using the purity formalism 
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to measure the polarization of various quark flavors us- 
ing identified hadrons in electroweak interactions, where 
parity violation provides polarization observables for free. 
We leave these possibilities to future work. We will fo- 
cus instead on one example of particular interest in spin 
physics: deep inelastic production of polarized hadrons 
(or jets). 

It is well known that polarized quarks can give rise to 
parity-odd correlations in the final state [ pi| . We illus- 
trate here how virtual photoproduction of vector polar- 
ized final states (see Ref. pq]) can be conveniently ana- 
lyzed using a formalism based on purities. Our results, 
written explicitly for virtual photon-nucleon scattering, 
generalize trivially for other hard hadronic processes. For 
dcfiniteness, we choose a helicity basis oriented along the 
virtual photon direction (+z) and neglect finite angle ef- 
fects, so that the fragmentation helicity basis is collinear 
with the quark helicity basis. We then combine the po- 
larizations p h z of chosen final states h (possibly specified 
by kinematics and choice of target) into a column vector 
P = {{p h })- O ur result is that at leading twist, for an 
unpolarized beam incident on a target with longitudinal 
polarization +z, 



where the matrix elements of AV are defined by 



AP 



(8) 



The AD^ are helicity difference fragmentation functions, 
which are assumed to be related, modulo dilution factors, 
to the quark helicity structure of the final state. Also, 
just as the purities satisfy the constraint Yli=f f = 1> 
the quantities AP/ 1 that enter Eq. (0) are normalized by 
the longitudinal polarizations p h z using a polarized beam 
(+i) and an unpolarized target: 

p h = AP/ 1 - (9) 

»=/,/ 

Equations (@)-(@) also hold for transversely polarized 
hadrons produced off a transversely polarized target, pro- 
vided we everywhere make the replacement A — > Aj_ 
(helicity ^transversity) . Indeed, formulas such as Eqs. 

are generic to processes with polarization observ- 
ables in the final state, where the (AD/D)^ are the corre- 
sponding analyzing powers. We therefore conclude that 
the purities previously defined are also the natural quan- 
tities for performing a global analysis of polarized hadron 
(or jet) production in hard processes. 
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